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We consider the problem of 'discrete-time persistence', which deals with the zero-crossings of a con- 
, tinuous stochastic process, X(T), measured at discrete times, T = nAT. For a Gaussian Stationary 

Process the persistence (no crossing) probability decays as exp(— 6dT) — [p(a)] n for large n, where 
a = exp(— AT/2), and the discrete persistence exponent, 6d, is given by 6d = Using the "Inde- 

pendent Interval Approximation", we show how 8d varies with AT for small AT and conclude that 
experimental measurements of persistence for smooth processes, such as diffusion, are less sensitive 
(— I I to the effects of discrete sampling than measurements of a randomly accelerated particle or random 

. walker. We extend the matrix method developed by us previously [Phys. Rev. E 64, 015151(R) 

(2001)] to determine p(a) for a two-dimensional random walk and the one-dimensional random ac- 
celeration problem. We also consider 'alternating persistence', which corresponds to a < 0, and 
calculate p(a) for this case. 

PACS numbers: 05.70.Ln, 05.40. +j, 02.50.-r, 81.10.Aj 
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I. INTRODUCTION 



Persistence of a continuous stochastic variable has recently been a subject of considerable interest among both 
theoreticians and experimentalists. Systems studied include randomly driven single degrees of freedom , simple 
diffusion from random initial conditions f|||], models of phase separation fluctuating interfaces [|l8p20|, and 

reaction-diffusion processes pl|-p4[. For a recent review see ^5|. Persistence is the probability, P(t), that a fluctuating 
non-equilibrium field, at a particular space point, has not crossed a certain threshold (usually its mean value) up to 
time t. In most systems studied, scale invariance implies that for large t, P exhibits a power-law decay, P{t) ~ t~ e , 
where the persistence exponent 9 is non-trivial due to the dependence of P(t) on the whole history of the system. 
Experiments have recently measured 9 for the coarsening dynamics of breath figures |2(| , liquid crystals |27j , soap 
bubbles |||] , and diffusion of Xe gas in one dimension |2^| . 

In this paper we consider the following problem: in any experimental (or numerical) measurement of 9, the stochastic 
variable studied, x(t), will have to be sampled discretely. It is therefore possible that x{t) could cross and then re-cross 
' its threshold between samplings, resulting in a false-positive classification of the persistence of x(t). If the sampling 
is logarithmically spaced in time (as was the case in p9|), then such undetected crossings will make the measured 
■ persistence exponent smaller than theory predicts (while if the sampling is uniform in real time only the prefactor 
is changed). This problem has been studied in |30| for the case of a random walker in one dimension, a simple 
Markovian, 'rough' (i.e. with a fractal distribution of crossings |3l)) process. Here we extend that work to consider a 
simple, non-Mar kovian, 'smooth' process: a randomly accelerated particle. Since most of the more complex processes 
studied experimentally are 'smooth', this paper is a step towards understanding the effect discrete sampling has on 
those measured persistence exponents. 
£j , There is yet another motivation for studying the persistence of a discrete sequence as opposed to that of a continuous 
process. It turns out that many continuous processes in nature are stationary under translations of time only by an 
integer multiple of a basic period (which can be chosen to be unity without any loss of generality). For example, 
the weather records have this property due to seasonal repetitions. It has recently been shown |]32f| that for a wide 
class of such processes, the continuous time persistence P(t) is the same as the persistence P(n) of the corresponding 
discrete sequence resulting from the measurement of the continuous data only at discrete integer points. In general, 
the calculation of the persistence of a discrete sequence is much harder than that of a continuous process, except in 
special cases where P(n) can be computed exactly [ [32|]33f| . The tools that have been developed over the last decade 
for studying the persistence of a continuous process are often not easily extendable to the case of discrete sequences 
and one needs to invoke new techniques, some of which are presented in this paper. 

The layout of this paper is as follows: In Section II we use the Independent Interval Approximation (IIA) |5|,|6| to 
find the first correction to 6 for small AT, where AT is the separation of samplings in logarithmic time (T = lnt). 
This gives us some indication of how significant the effect of discrete sampling is. For example, for the random walk 
9d — 9 cx — VAT and so the discrete exponent, 9d, and its continuum limit, 9, begin to deviate markedly as soon as 
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AT > 0. For the random acceleration problem, one finds 6n — 9 oc —AT, a weaker dependence on AT, while for the 
diffusion equation from random initial conditions, dp — 9 oc — (AT) 2 , so the effect of discrete sampling for small AT is 
much less significant in this case. These three systems display an increasing 'smoothness' of the underlying process, 
a concept we will expand upon below. 

In Section III we illustrate our general method for solving low-dimensional discrete persistence problems by consid- 
ering the case of a random walker in one-dimension, x(t) = rj(t), where rj(t) is Gaussian white noise. A brief account 
of this work was given in We map the problem to a Gaussian Stationary Process (GSP) in the variable X — xj\fi 
by transforming to logarithmic time. In the new variables, the random walk is represented by an Ornstein-Uhlenbeck 
process. We re-pose this problem in terms of a backward Fokker-Plank Equation (BFPE) and give the solution for 
the continuum case where the persistence exponent has the well known value 9 = 1/2. We then formulate the discrete 
persistence problem in terms of an eigenvalue integral equation. Employing a power series expansion of the integrand 
reduces the problem to a matrix eigenvalue equation, whose largest eigenvalue gives us the persistence exponent 
measured by discrete sampling. The concept of alternating persistence, in which the consecutive measured values of 
X lie on alternate sides of the threshold value, is also introduced, and 8d calculated for this case. 

Having demonstrated the matrix method on a single variable problem, we then apply it to two two- variable problems: 
the two-dimensional random walk, x(i) = r)(t), in a wedge geometry (Section IV), and a randomly accelerated particle 
in one dimension, x(t) = T)(t), which is a simple example of a non-Markovian process (Section V). Surprisingly, for 
alternating persistence in the continuum limit (AT = 0) the asymptotic probability of surviving one further sampling, 
p, is non-zero, a phenomenon contrary to our earlier suggestion [ p0| that this should only occur for 'rough' processes. 
Using scaled variables and logarithmic time, the random acceleration problem maps onto a damped simple harmonic 
oscillator, which we study for the overdamped case using the matrix method, finding 9jj as a function of AT. Using 
the results of a correlator expansion developed in [p4| we also study the underdamped case, an example for which the 
correlator is oscillatory (Section VI). We show that when the time interval, AT, between measurements is equal to 
the period of the oscillation, the problem is identical to the Ornstein-Uhlenbeck process studied in Section III, while 
for AT equal to one half of a period the problem reduces to that of alternating persistence in the Ornstein-Uhlenbeck 
process. The paper concludes with a summary of the results. 



II. SMALL AT CORRECTION TO 6 

The Independent Interval Approximation (IIA) ]||| uses the assumption that the intervals between zero-crossings 
of a GSP are independently distributed. Although this assumption is not valid for most processes, it nevertheless 
gives remarkably accurate estimates for 9 in many cases. Here we use the same assumption to find the first correction 
to 9 due to discrete sampling with a spacing in logarithmic time of AT for AT small. 

As AT is increased from zero, the first correction to 8 comes from paths which are always positive apart from 
one undetected double crossing between consecutive sample times. Let the probability of one such double crossing 
occurring in the interval T — nAT be PdUe(n, AT). Then, for T large, the probability, Pjj(T), that the stochastic 
variable is positive at all n samplings is given, to lowest order in AT, by 

Pd(T) = P a (T) + P dble (n, AT) (1) 

where Po{T) is the continuous-time persistence probability, given by Pq(T) ~ e~ 9T for T large, and Pdbie(n, AT) is 
given by 

P dble (n, AT) oc n e' 9nAT / dT x / dT 2 P X {T 2 - T x ) (2) 

JO JTx 

where Pi(T) is the probability distribution of the interval size, and we have assumed that the durations of different 
intervals are statistically independent. The latter assumption is precisely the IIA. It is clear from Eq. (^) that, to 
leading order in AT, we only require the form of P\{T) in the limit T — > 0. The function P\(T) can be found using 
the IIA and for the processes currently under consideration - the random walk, random acceleration and diffusion 
from random initial conditions - the small-T results are Pi(T) oc 1/vT, 1, and T respectively. All three cases are 
incorporated in the general form Pi(T) oc T a , with a — —1/2, and 1 respectively. Using this form in Eq. (|^) gives 

P D (T) = Po(T)+ in £T a+2 e- eT (3) 

where 7 is some constant. Since T — tiAT and Po(T) ~ e~~ eT , we have 

P D (T) = Ae- eT (l + BT£T a+1 ) (4) 
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where A, B are constants, and so, to lowest order in AT, 

P D (T)^Ae- eT e BT ^ a+ \ (5) 

Since P D (T) ~ e _eDT , we obtain 

9 D =9-B/XT a+1 . (6) 

For the random walk, random acceleration and diffusion from random initial conditions, 9d ~ 9 oc — V AT, —AT, and 
— (AT) 2 respectively. From this we expect that the discrete sampling is important as soon as AT is non-zero for the 
random walk, this being related to the fact that P\{T) — > oo for T — * 0, i.e. that the distribution of crossings is 
fractal and hence this process is 'rough' |n| . The persistence exponent for the random acceleration process is linear 
in AT and thus we expect it to be less affected by the discrete sampling for small AT . In sections III and IV these 
expectations will be confirmed using the matrix method perturbative expansion in a = e ^^/ 2 about a = 0. Finally, 
for diffusion from random initial conditions, the (AT) 2 dependence indicates that discrete sampling will be relatively 
unimportant for small AT. Finally we note that, although Eq. (^) has been derived using the IIA, we expect it to be 
valid quite generally. In particular, the a = 2 result for P\{T) has been proved correct by Zeitak |37), the IIA even 
giving the correct coefficient. 



III. RANDOM WALK IN ONE DIMENSION 



Let us consider the simple case of a random walker in one dimension, 

±(t) = V (t) (7) 

where r](t) is Gaussian white noise with zero mean and (r](t)r)(t')) — 2D5(t — t'). For convenience we map this process, 
which is non-stationary in time, to a Gaussian Stationary Process (GSP) by defining new space and time variables 
X = x/Vt, T = Int. Eq. (g) then reads 

^ = _ M X(T) + r ? (T), (8) 

where (r)(T)r)(T')) = 2D8(T — T") and fi = 1/2 for the random walker, although other values of /j, can be considered. 
This Ornstein-Uhlenbeck problem can be more usefully solved using the BFPE. Let Q(X,T) be the probability that 
the random walker, starting from X at time T = 0, does not cross the origin (X — 0) up to time T. Then from (||) it 
can be shown that Q(X,T) satisfies 

dQ/dT = D d 2 Q/dX 2 - fiXdQ/dX , (9) 

with initial condition Q(X, 0) = 1 for all X > and boundary conditions Q(0, T) = and Q(oo, T) = 1. The solution 
is 



Q(X,T) = Erf 



-.X 



(10) 



where Erf(x) is the error function and D 1 = D/fi. For large T (and positive fi), P(T) <~ Ae _AlT , which corresponds 
to P(t) ~ t~ e in real time, with 9 — /i. 

We now consider discrete persistence, i.e. the probability Q n (X) that starting at X our field/ variable is positive 
at all the discrete sample times T\ = AT, — 2 AT, T n = nAT. This is relevant for experimental or numerical 
determinations of 9 since in practice one will have to sample only at discrete points. In ref. p9| , for example, the 
sampling is done logarithmically in real time (which corresponds to sampling uniformly in logarithmic time). Note 
that 0d will in general be smaller than the continuum value 9 since any even number of crossings between samplings 
will go unnoticed. One can write down a recurrence relation for Q n (X): 

POO 

Q n+1 (X) = / dYQ n (Y)P(Y,XT\X,0) (11) 
Jo 

where P(Y, AT|A, 0) is the Greens function, i.e. the probability that a particle starting at X at time will be at Y 
at time AT. For a Gaussian process, P(Y, AT|X, 0) can be found from the mean and variance of X(T), 
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P(Y, AT|X,0) 



exp 



(Y-aX 



,2 1 



2D'(1 - a 2 ) 



where a = exp(— /uAT). This gives us the discrete analogue of the BFPE <ffi), 



WW 



dFQ„(y) 



^vrL^l - a 2 ) 



exp 



2L>'(1 - a 2 ) 



(12) 



(13) 



with the continuum equation being recovered in the limit AT — > 0. 

Making the change of variables a; = X/ yjD'{l — a 2 ), y = Y j y/D'{l — a 2 ), and Q ra (-^) = <2n( x ) gives, 

1 f 00 

Q'n+x^) = j= J Q dyQ' n (y)cx P {-(y~ax) 2 /2} (14) 

At late times, with fj, > 0, we expect Q' n+ i{x) — pQ' n {x) in analogy with the continuous case where P(T + AT) 
P(T)e- eAT , so p = e- SDAr . We therefoi 
an eigenvalue integral equation for q(x) 



P{T)e eAT , so p = e eoAr . We therefore expect that, for large n, Q n (x) — ► p n q{x). Substituting this into ( |l4| ) gives 



2tt Jo 



dyq(y)evp[-(y - axf /2] , 



(15) 



with an eigenvalue p(a) that depends continuously on a. Eq. ( |l5| ) has an infinite number of eigenvalues, but at late 
times only the largest will remain since J2i Pi ~ Pmax f° r l ar g e n (there is no other eigenvalue contiguous to the 
largest eigenvalue). By symmetrising the kernel in (|l5|), one can use the variational method to find a rigorous lower 
bound for p ]30[ |. This method cannot be applied to the random acceleration problem, however, since the kernel 
cannot be symmetrized. For AT — > oo, correlations between the nth and (n+2)th samplings are negligible, so 



P n w [Prob (two consecutive points have the same sign)]' 
1 +i(sign[X(0)]sign[X(Ar)]) 



= (i + ~~arcsin[C(AT)] 



(16) 
(17) 

(18) 



where CiT^ — T\) = exp[— p,(Ti — T\)] is the normalized autocorrelation function of X(T), i.e. C(T) = 
{X(T)X(0))/(X 2 ). Hence, for large AT, 



1 a 

P= 2 + ^ 



Eq. ( |l5| ) can also be solved perturbatively by expanding the exp(axy) term as a power series in a, 



pq(x) 



dyq(y)e 



-y 2 /2 p ~a 2 x 2 /2 



E 



(axy) r 



Defining 



gives 



2/2 poo 



dyq(y)y m e-y 2 / 2 



(19) 



(20) 



(21) 



pq{x) 



-a 2 x 2 /2 



2tt 



E 

m=0 



2 /o 

Multiplying through by e~ x 1 x n and integrating over x > gives 

oo 

pb n — ^ ^ Anm 



m=0 



(22) 



(23) 
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where 



At. 



^4^(1 + a 2 ) 



2 a 



1 



(n+m)/2 



r[(n + m + l)/2] 
V 'n\m\ 



(24) 



By this method, the eigenvalue integral equation (|15|) has been converted to the eigenvalue matrix equation fl23|), 
and the problem reduces to computing the largest eigenvalue of an N x N submatrix whose (n,m)th elements decrease 
exponentially in n + m. In pcj we determined p numerically to one part in 10 12 . We have also algebraically found 
Pmax as a series expansion in a, the first four terms being: 



a 

IT 



48 - 36tt + 7tt 2 

67T 3 



a 3 + 



(25) 



The coefficients up to order a 49 are given in appendix 2. For AT — > 0, a — > I and convergence becomes progressively 
slower. However, the variational method still works in this region. 

One may also consider the case of alternating persistence, i.e. the probability that X(n/KT) is positive for every 
even n and negative for every odd n (or vice versa). The limit of integration in Eq. ( [l5| ) then changes, giving 



pq{x) 



1 



V2^ J- 



dyq(y)exp[-(y - axf /2] 



(26) 



Substituting y — > — y, swapping the limits of integration and using q{—y) = q{y) (since the process is symmetric 
around y = 0) gives 



pq{x) 



dyq(y)exp[-(y + ax) 2 /2] 



(27) 



which is identical to (|l5|), but with a replaced by —a. Therefore replacing a by —a in the matrix equation (|2^) will give 
the alternating persistence eigenvalue p a it- The case \a\ > 1 may also be considered, this corresponding to p < and 
hence to an unstable potential in the Ornstein-Uhlenbeck process. For the alternating case, a < — 1, the calculation 
proceeds as before. In fact, from Eq. ( p4| ) it can be seen that p(l/a) = p(a)|a| for a < 1, so one need only investigate 
alternating persistence in the range — 1 < a < 0. For a > 1 it is possible for the walker to escape to infinity, i.e. 
qn(x) y^> for n — > oo. The asymptotic limit for q(x) when a > 1 may be found using the matrix method in the 
following way. Since q(x) — » 1 as x — > oo, it is more convenient to study a new function u{x) defined by the relation, 



q{x) = 1 - q(0) / u(y)dy 



(28) 



where q(0) is fixed by Eq. ( |28| ) with x — 0. Substituting q(x) from Eq. (f28|) into Eq. ( |15| ) (with p = 1 since we are 
finding the stationary state) we find, after some algebra, that u(x) satisfies the integral equation 



u{x) 



'2tt 



u(y)e 



-(y-ax) 2 /2 



dy 



(29) 



Note that, unlike Eq. (|15|), which determines q{x) only up to an overall multiplicative constant, Eq. (|29J) is an 
inhomogeneous equation which fixes u{x) absolutely, as one expects on physical grounds. 
As before, we expand the factor exp(axy) in ( |29| ) as a power series to obtain, 



u(x) 



-a 2 x 2 /2 



2n 



dyy n u(y)e-y 2 / 2 



(30) 



Multiplying through by x m a ^ m ^ e s integrating over positive a;, and defining 

ft n/2 



71! Jo 



dyy n u(y)e-y ' 2 , 



(31) 



gives 
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c n — d[A n Q -\- ^ ^ -A nm c m ] , 



(32) 



where A nm is our previous matrix given by Eq. (p4), and 



u(x) 



Eq. (B2) can be solved by matrix inversion: 



n=0 



(33) 



(34) 



where 



(35) 



The solution converges rapidly as a function of the size, N, of the matrix. In practice, N of order a few hundred gives 
very precise results. 

Figure |l| shows p(a) for both alternating and normal persistence, while Figure ^| shows 0(a). Note that even for 
a = 0.96, Ojj is significantly below the continuum value of 1/2. We argued in section II that the difference 8(1) — 9(a ) 
decreases at AT 1 / 2 , i.e. as (1 — a) 1 / 2 for a — > 1, implying a square-root cusp at a = 1 in Figure ^. Wong et.al. P9| , 
measuring persistence in one-dimensional diffusion of Xe gas, sampled logarithmically in time such that in log-time 
their AT was about 0.24. If the process were a random walk, this would give a ~ 0.9 and a difference between 9o 
and the continuum 9 of about 20%. For the diffusion equation that describes the experiment, however, the approach 
to the continuum is more rapid (see the discussion in section II), and rather accurate results arc obtained even for 
a = 0.24 H§. 

In Figure [| we show the eigenfunction q(x) for the case a — 0.5. This function was obtained by substituting the 
eigenvector of the matrix A corresponding to the largest eigenvalue into Eq. (|2^). The asymptotic behavior for large 
x is q(x) ~ x u , where v = lnp/lna ~ 0.530661 for a = 0.5 [|o| . Figure || also contains a plot of q(x), given by Eq. 
(pq), for the case a = 2, which corresponds to an unstable potential. 



IV. RANDOM WALK IN 2-DIMENSIONS: WEDGE GEOMETRY 



Having illustrated the perturbative method on a simple Markovian case for which another approach (the variational 
method) is available, we intend to study a simple 'smooth' non-Markovian process, the random acceleration problem, 
x = r/(t). This process is equivalent to a Markovian problem in two variables, x(t) and v(t), where v = 77(f) and x = v. 
Before dealing with this problem, we will first consider another two- variable Markov process, namely the random walk 
in two dimensions, x — rj x (t), y = rj y (t), with (rji(t)) = and (r)i(t)rjj (t 1 )) = 2D5ij5(t — t'). Using the perturbative 
approach on this pedagogical problem will clarify its use on the random acceleration problem. 

Consider a wedge of angle a whose boundaries are absorbing, let our random walker start inside this wedge at radial 
position r and angle <f>, with < <p < a. Making the change of variable R = rj\ft, T = kit, converts tke problem into 
a GSP, as in section III. The corresponding BFPE, i.e. the two-dimensional analog of Eq. (||), is 

dQ/dT = V 2 Q- pRdQ/dR (36) 

where ji = 1/2 for the random walk, though we will treat /i as arbitrary, and Q(R, 4>, T) is the survival probability of 
the particle at time T given that it started at (R, <fi). The initial condition is Q(R, <t>, 0) = 1 for R > and < </> < a, 
while the boundary conditions are Q(R, 0, T) = 0, Q(R, a, T) = 0, Q(0, 0, T) = 0, and Q(oo, 0, T) = 1 for < cj) < a. 
The solution can be obtained using separation of variables. Its asymptotic form for T — > oo at fixed R and cj) is 

Q(R, <p, T) oc R 7 */* sm(-K(j)/a) exp(-^7rT/a) , (37) 

giving 9 = p/Kja. Now consider the discrete persistence. The analog of equation (^) is, 

pq(r) = ^f dv'q(v') exp[-(r' - ar) 2 /2] , (38) 
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where a = exp(— /xAT) as before, r = (1 — a 2 ) 1 / 2 R, and the integration is over the wedge. Note that for a = n one 
can do the integration over x and recover the one-dimensional result. In polar coordinates, Eq. ( |38|) becomes: 



pq(r, 4>) = — I r'dr' I d(f>'q(r', ft) exp[-(r' 2 + aV - 2arr' cos(0 - ft))/2] • (39) 







2^ Jo 

As before, we expand the exponential term containing the mixed terms, 

exp[arr' cos(</> -ft)] = am ' + ™' (rr') m ' +n ' [cos{ft) cos(0')P' [sm(0) sm(ft)] n ' /[m'ln'l] , (40) 



and define 



(m'+n')/2 /,/37j3j 



a />oo 



'/ drV m ' +,i ' +1 e~ r ' / 2 [cos m '(0')sm™'(0')]9(^,0') , (41) 
o Jo 



giving 

p?(r, 0) = [e- a2r2/2 /(2ir)} ^ r m ' +n 'a (m ' +n ' )/2 [cos m '(0) sin"' fc/ViiiW . (42) 
Multiplying through by 

L(™+»)/2/V^!^ r m+«+l e -r 2 /2 cos ™(^) sm "(0) , (43) 

and integrating over the wedge, gives 

P^ran — ^ ^ ^mn,m' 'n'^m' ' n 1 • (44) 

Note that, by regarding mn as a single index, we can treat A mn ^ m i n i as a conventional matrix when solving numerically 
for p. The matrix elements are, 



2a 



-| (m+n+m'+n' +2)/2 

r[(m + n + m' + n' + 2)/2]J m „, mW (45) 



4ira\/ mlnlm'ln'l 
where 

W»'= / a #[cos(^] m + m '[sin(0)]" + "' . (46) 
Jo 

The problem has been reduced to an eigenvalue matrix equation in four variables which we solve numerically. We 
have also found p as a series expansion in a. As before we can consider alternating persistence, i.e. the probability 
that between each sampling the particle alternates between < 4> < a an d tt < 4> < ot + ir. Figure || shows a plot of 
the results. The case a = n corresponds to the one-dimensional case while a = and a = it trivially give p=0 and 
p=l respectively. 

For AT — ■+ oo, p — > a/2ir since the system re-equilibrates between samplings: In this limit the function p(a) is a 
straight line from (0,0) to (27r,l). 

In the limit a — ► 0, we can expand the ft integral in ( |39| ) to first order in a, and set <fi = = ft to obtain 

o If 00 

-q(r,0) = — / r'dr'q(r',0)exp[-(r' 2 +a 2 r 2 -2arr')/2] , (47) 
a 2tt Jo ' 

which is identical to the one-dimensional random walk equation (Jl^) but with r'dr 1 rather than just dr' , and p replaced 
by y/2~Trp/a. This gives a matrix equation similar to that for the d = 1 case, but with p oc a. Also p is a monotonically 
increasing function of a as in the d = 1 random walk case. The same argument applies to the alternating case (where 
a — > — a). This explains the linear dependence of p on a for small a in Figure |J and the qualitative a-dependence of 
the slope at a = 0. 

In the limit a — *■ 2-7T, we expect p — ► a/2-n since the "mean free path" of the particle between samplings is very much 
greater than the width of the absorbing region, so that the particle does not "see" the absorbing region and hence the 
probability distribution is not affected by it . Thus doubling (2tt — a) doubles the probability of being absorbed. This 
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is true for any AT although since the "mean free path" cx V AT, as a is decreased from 2tt the linear relationship will 
fail sooner for smaller AT. The same argument applies for the alternating case, but since there are two alternating 
absorbing regions, the probability distribution has more time to equilibrate near each absorbing region, thus the linear 
dependence breaks down at smaller a for a given AT. This explains the linear dependence near a — 2ir in Figure ^, 
and its breakdown as a decreases. 

In the continuum limit (a = 1) we have the result 9 = (iir/a. Since p = a e ^, for a = 1 we get p = 1 for all a. The 
a > plots in Figure || are tending to this limit for a— ► 1. In the alternating continuum limit (a = —1), the particle 
cannot survive for a < tt, since the two sectors that the particle must occupy alternately are disjoint. For a > tt 
these two sectors overlap, and it is clear that the alternating persistence problem (a = —1) is equivalent to the usual 
persistence (a = 1) but with a replaced by a — tt, i.e. = irp,/ (a — tt). So p = for a < w, and p = 1 for a > tt, i.e. 
we get a step function at a = tt. The a < plots are tending to this limit for a — > — 1. 

From equation flt5|) it can be seen that p(l/a) = a 2 p(a) for a < 1, which enables one to construct the results for a 
in the range a < — 1 from those in the range — 1 < a < 0. The extra factor oil/ a compared to the d — 1 case is due 
to the phase space factor rdr rather than dx in the eigenvalue integral equation. 



V. RANDOMLY ACCELERATED PARTICLE 



We now consider the non-Markovian random acceleration problem. We first recast the equation x = rj(t) as the 
two first-order equations v — rj(t) and x = v. After the change of variables V = v/t 3 / 2 , X = x/t 1 / 2 , T = hit these 
equations become 



dV/dT = -aV + v(T) 
dX/dT = -/3X + V 



(48) 
(49) 



with a — 1/2, p = 3/2 for the random acceleration problem, the process being 'smooth' for any a, p. The noise 
correlator is (rj(T)Tj(T')) = 2D5(T-T'). By eliminating V, Eqs. @, @ can be written as 



d 2 X/dT 2 + (a + P)dX/dT + afiX = -q(T) , 



(50) 



which is the e qua tion for an overdamped simple-harmonic oscillator, a — P corresponding to critical damping. Note 
that equation ( |50| ) is symmetric under a <-> p. Letting T — > T/(a + p) and X — + Xj (a + P) 3 ^ 2 gives, 



X + X 



a/P 



(1 + a/py 



-X = V (T) 



(51) 



showing that 9 = af(a/p) and so only differing ratios of a to P need be considered. 

For continuum persistence, Sinai pj| and also Burkhardt 0] have shown that = 1/4 for the random acceleration 
problem (a — 1/2, P — 3/2). For ratios other than P = 3a, no analytic solutions have been found, except for the case 
a/j3<l which we give in this paper. 

To derive our usual eigenvalue integral equation, we need to find the propagator P(X, V, T + AT|Y, U, T), that is, 
the probability of going from (Y,U) at time T to (X,V) at time T + AT. For a Gaussian process, finding the means 
and variances will completely specify the distribution. 



Let X = X- < X > , V = V- < V > and define A = e 

1 



B = e /3AT '. Then the propagator is, 



P(X,V\Y,U) 



where 



27iVDetA/ 



M 



exp 



±{X,V)M-\X,V) 



(52) 




IxvX 



and 



X = X-YB 
V = V -UA 



u 



P-a 



(A-B) 



(53) 
(54) 
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(* 2 ) = a pm_a)Vp + a} ^ 13 ~ ^ ~ ^ + ^ + 4aPAB _ a(/3 + a)B2) (55) 
(xv) = ^-l— W - a) - 03 + a)A' + 2aAS) (56) 

l/ 2 \ = -(1- A 2 ) (57) 



a 

and the eigenvalue-integral equation is, 



pP ao (X,V) = dY dUPUY.^PiX.V^U) , (58) 

JO J-oo 

where we are working in the forward variable P ao (X 1 V) 1 the probability of finding the particle at position X with 
velocity V given that it has been found at positive X at all previous samplings. Making the change of variable, 

= x ' 73537 = v < 7353 = V ' 7353 = u > P -( X ' F ) = ^ w ) S ives: 

pf(x,v) = -VDctM / dy / d w /(y, u ) e -( £2 < v2 >- 2i5 < xv >+" 2 <* 2 >) (59) 

JO J-oo 

Substituting 

v = v — Au, 
A — B 

x = x — yB — u— , (60) 

p — a 

into Eq. (|59|), the exponent will contain all pair combinations of x,v,y,u. This kernel is not symmetrisable in x <-> y 
with v <-> u, hence the variational method is not applicable. Expanding in the four mixed terms, xy,xu,vy,vu, and 
using the same method as before, we obtain a matrix equation of the form 

pl c<d = M c> d, e jl e ,f (61) 

(see appendix 1). As previously, we have found p numerically and also as a series expansion in a = exp(— AT/2). The 
results for p(a) are presented in Figure ||, and the corresponding results for 9(a) in Figure ||. The coefficients of the 
expansion up to order 25 are given in Appendix 2 for the case a = 1/2, f3 = 3/2. As for the random walk case, the 
series have not yet converged for a — * 1. Since we now have a four-dimensional matrix, the problem is more apparent 
since we have only been able to reach 0(a 25 ). We extend the curves by use of Pade Approximants p^j , which relies on 
p(a) being 'smooth', and so do not work for the random walk for a — ► 1 because of the \J AT cusp. Also, (see Figure 

, the Pade method does not work in the alternating case for large (3 because of the sharp downturn for a — * — 1 . 
Nevertheless, in the remaining cases it significantly extends the valid range for p, and we use this for plotting the 
remaining alternating persistence cases. However, since 9 — iln(p)/ln(a), when a -> 1 the l/ln(a) accentuates the 
slight error in p. To remedy this, we add an extra term to the Pade polynomial in the numerator (or denominator) to 
enforce the constraint p(l) = 1. Figure || shows plots of the eigenvalue against exp(— AT/2) for a = 1/2 and various 
values of (3. Also plotted are the alternating persistence results. 

Surprisingly, for the continuum limit of alternating persistence p -/-> 0. Indeed, p m 0.108 for all a, (3. Previously, 
we had thought [j30[ that non-zero p could only occur for a 'rough' process (which has an infinite number of crossings 
in finite time). The reason for these results is as follows. For alternating persistence, X (and similarly V) must cross 



zero inside a time AT so for AT — > 0, X, V must both tend to zero. As a result, Eqs. (Eq), (49) become V(T) = rj(T), 
X(T) — V(T), or equivalently X — r](T), thus removing the a, (3 dependence. This latter equation is invariant under 
changes in T and hence AT, provided we rescale X . It therefore gives a non-zero p even for AT — > 0. This implies that 
any process whose equation is (or becomes, for AT — > 0) time-scale invariant has p a it ^ for AT — > 0. For example, 
the diffusion equation from random initial conditions is equivalent to the n — > oo limit of the process d n x/dt n = rj(t) 
and will reduce to d n X/dT n = n(T) for AT — > 0. 

We can compare our results to exact results for two limiting cases. For AT — ► oo, as before we expect p = 
5 + — arcsin[C(AT)], where for this case the normalized correlation function is 

C(AT) = El ™ (62) 

p — a 



9 



Hence p =1/2+ i^e~ oA/ + 0{e~ 2a ^ , e~P m ) for /3 > a. For a//3 -» 0, the normalized correlator, Eq. (|6g), 
becomes C(AT) = exp(— aAT), the random walk correlator. In fact, Eq. (|50| ) reduces to the random walk equation, 
X = —aX + rj(T), if the limit a <§C (3 is taken after the change of variables X = \/f3\/a~, T — r/a, when it becomes 
clear that the inertial term is negligible for a/ [3 — > 0. For a//3 <C 1, we can also solve for in the continuum limit, 



AT = 0, by expanding around the Markov process X = -aX + ri{T) to first order in y/a/P HHJlJ. The correlator 
C(T) can be written as 



C{T) = e- al + 



(3-a 



-oT _ „-/3T 



(63) 



Using the standard perturbation expansion 1 36 1 for for a process with correlator C(T) = e~ aT + ea(T), with e < 1, 



a 1 



gives 



2a a 



a 1 - 



a(T) 



("J _ e -2aT^3/2 



+ 0(e 2 ), 



/7F /3 — a 



2a 



/r 



/3- 



2a 



To first order in y/a//3, therefore, 



= a 1- 




(64) 



(65) 



(66) 



For a = 1/2, = 24, this gives 9 = 0.442 + 0(a/(3), which is consistent with the value of = 0.429 obtained from 
the constrained Pade. 

The series for p(a) in powers of a agree with those found using the correlator expansion |34| for the various values 
of a and j3 tried. This was used as a powerful check of the accuracy of the correlator expansion. 

As before, the eigenfunction P ao (X, V) can be reconstructed from the eigenvector corresponding to the largest 
eigenvalue using equation (A4). The result is displayed in Figure [7] for the standard case a — 1/2, (3 = 3/2. 



VI. UNDERDAMPED NOISY SIMPLE HARMONIC OSCILLATOR 



We now consider the underdamped case of the noisy simple harmonic oscillator. The persistence exponent for this 
case is not known even in the continuum limit. Furthermore the correlator of the process is oscillatory, which to our 
knowledge has not yet been studied. From equation ( |5C| ) it can be seen that the complex values 

a = 7 + ilu 

/3 = j-iu (67) 

correspond to an underdamped oscillator. The corresponding equation of motion is 

X + 2 7 A + ( 7 2 + u?)X = n(T) (68) 

where the angular frequency of oscillation is lo and the decay rate is 7. Substituting a and (3 into the correlator, Eq. 
@, gives: 

C(T) = exp(— jT) [cos(wT) + - sin(wT)] . (69) 

Since a, [3 are complex, we choose to study this process using the correlator expansion about /XT — > 00 which was 
developed in Q. We merely substitute the correlator, Eq. (|69|), into our 14th order series expansion for p{/XT) and 
hence find 9d- Note that for the random walk and random acceleration problems above, 14th order corresponds to 
order a 14 . As before, 0d = 7/(^/7) so we choose to keep 7 = 1/2 and vary u>. Figure || shows plots of 9 against 
XT for various values of u>. Also shown are the random walk and persistence and alternating persistence exponents, 
which are numerically identical to the underdamped SHO when luXT — 2mn and luXT — (2m + 1)tt respectively, with 
m an integer. Note that, as before, the series have not converged sufficiently for small XT. 
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The (at first sight surprising) identity of the random walk and underdamped SHO persistence exponents for certain 
values of AT is interesting, and the explanation rather simple: Since these are Gaussian processes of zero mean, their 
correlators possess all the information about them. When studying discrete-sampling persistence, we are in effect 
studying a sequence of random variables, X\, X2, X3, ...A„, whose correlator we know since (X p X q ) = C(\p — q\ XT). 
Any processes which have identical correlators for all C(mXT), m integer, will have identical discrete persistence 
exponents for this value of XT. From equation (|69| ) it can be seen that for loXT — 2rmr the correlator reduces to 
the random walk correlator, exp(— 7 AT), hence the agreement of the persistence exponents. For ojXT = (2m + ljir, 
C(mXT) = exp(—"fXr)(—l) m , that is, the correlator alternates in sign. This corresponds to alternating persistence 
since X(mXT) has in effect changed sign for m odd. 

At this point we make a general comment about the variation of discrete persistence exponents with XT. Consider 
a process for which we know 9(XT) for some range of AT, e.g. from using the matrix method or correlator expansion 
for small AT: If we change XT to XT/m, m integer, we will sample at exactly the same times as before, plus the 
intermediate times. Hence some paths that were persistent before will be excluded, and so the asymptotic survival 
probability decreases. Thus, 

e D (XT/m) > MAT). (70) 

This of course does not imply that #d(AT) is monotonic. Taking the limit m — > 00, and assuming that 9d(XT) is 
smooth for AT — > 0, we get the result that 

9>9 D (XT), VAT. (71) 

This provides a lower bound on 9 for when our large XT expansions do not converge up to the continuum limit. We 
can use this to see that, for 7 = 1/2, 9 > 1.22, 2.60, 4.49, and 7.19 for uj = 2, 4, 8, and 16 respectively. 



VII. CONCLUSION 



In this paper we have extended our earlier treatment of the discrete persistence exponent for the random walk 
with an absorbing boundary at the origin f30| to the two-dimensional random walk with absorbing boundaries on a 
wedge, and the random acceleration process with an absorbing boundary at the origin. While the latter is perhaps 
the simplest continuous-time non-Markovian process, both of the processes discussed in this work have the simplifying 
feature that they can be written as a Markov process in two dimensions. We have shown that, in both processes, 
the discrete persistence probability after n measurements, Q n (x), for a process starting at x, has the asymptotic 
form Q n (x) ~ p n q(x), where p and q(x) are the largest eigenvalue and corresponding eigenfunction of a certain 
eigenvalue integral equation. They have been evaluated to high precision by converting the integral equation into a 
matrix eigenvalue problem, from which high-order series expansions in powers of a — exp(— pXT) have been obtained, 
where AT is a uniform measurement interval for a Gaussian stationary process (GSP). The random walk and random 
acceleration problems have been mapped onto GSPs in logarithmic time, T = lni, and for these (and similar) processes 
the calculations presented here apply to the case of measurements uniformly spaced in T. 

The case of alternating persistence, in which the measured values of the stochastic variable take positive and 
negative values alternately, has been discussed using a formal continuation of a to negative values. The case a > 1 (or 
a < — 1 for alternating persistence), corresponds, for a GSP, to motion in an unstable potential. For a > 1 there is a 
non-zero probability, q(x), that the process, starting at x > 0, is never measured to be negative, and we have shown 
how to calculate it, with explicit results (see, e.g., Figure 3) for the Ornstein-Uhlenbeck process. 

For the random acceleration process, the corresponding GSP is a noisy, overdamped harmonic oscillator, X + (a + 
(3)X + a(3X = r](T), with a = 1/2 and (3 — 3/2. This process is clearly of interest for other values of a and (3, since 
even the value of 9 in the continuum limit (AT — > 0) is not know exactly except for (3 = 3a. We have obtained a 
perturbative result for the continuum limit in the limit a -C (3. We have also investigated the underdamped case, 
corresponding to complex a and (3, and shown that for discrete measurements with time step AT equal to an integer 
number of oscillation periods, the persistence properties are identical to those of a random walk, while for an odd 
half-integer number one recovers the alternating persistence of a random walk. 

The methods presented here become increasingly unwieldy as the order of the stochastic differential equation 
increases. Recently, a power series approach has been developed in which the eigenvalue p is expanded in powers of 
the correlator, C(fcAT), evaluated at integer multiples, k, of the time step between measurements, with the maximum 
value of k depending on the order of the expansion 34 1. While not as powerful as the matrix method for systems 
described by low-order stochastic differential equations, the series approach has the advantage that it can be applied to 
any GSP, including those (such as diffusion from random initial conditions) which cannot be described by a differential 
equation of finite order. 
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The 'simple' persistence problem discussed here deals with the probability of detecting no zero crossing in n measure- 
ments. The statistics of the number of zero crossings, i.e. the probability to observe m crossings in n measurements, 
is also of interest. Results obtained by applying both the methods of this paper and the series expansion approach of 
ref. (34] will be reported in a separate publication. 

This work was supported by EPSRC (UK). 
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APPENDIX A: CALCULATION OF RANDOM-ACCELERATION MATRIX EQUATION 

Let the xv coefficient in the exponent be a xv etc, then: 
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2(xf) 



-B 2 (v 2 

A-B I - 



2B 



(V 2 

f3-a \ 



2AB ljtv\ 



-A 2 (X 2 
2B(v 



A-B\ /~o\ nA A- B 
(V 2 ) + 2 A . 



(xv 



-2B (xv\ 



Expanding the exponent of Eq. ( [59|) in powers of terms which mix (x,v) with (y,u), we get, 

2 poo />oo 00 

pf(x,v) — -VDetM / dy duf(y, u) — {a xy xy + a vy vy + a xu xu + a vu vu) r 
n Jo n =o n - 



x exp[-(a xx x 2 + a xv xv + a vv v 2 )] 



Letting 



gives 



00 poo 



I c ,d= J dx I dv f(x,v)x c v d exp[-(a xx x z + a xv xv + a vv v 2 )} 



pf(x,v) 



7r * — ' ' — ' fe+ f)\ \ e I * — ' ' — ' \r I \ s 

e=0 /=0 v ' v 7 r=0 s=0 



e\ (f 



a e ~ r a r a f ~ s a s x e ~ r+f ~ s v r+s 

u X y u vy u ' X u u mi^ u 



x exp[-(a xx x 2 + a xv xv + a vv v 2 )] 
Multiplying through by x c v d , and integrating over x > and over all v, gives, 



pled = - VDetM G c , d ,e,fl, 



where 



and 



giving 



Gc ' d ' eJ ~ ( e + /)! ( e ^ ^ (0 a *» r °"i' * u Sa S«-°e-r+/- s +c,r+ s +d 



D a>b = dx dvx a v b exp[-(Av 2 -Bxv + Cx 2 )} , 



^ = \/Ie 
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where [|] indicates the integer part of | and 

B = a xv + a yu (A10) 

C = G^a; Gyy ■ (All) 

Thus the problem has been reduced to computing the largest eigenvalue ofaniVxTVxiVxTV operator. As in the 
wedge case, G c ,d,ej decouples into G( c ,d),( e ,/)' For alternating persistence, the range of integration over y in Eq. ( [A2] ) 
should be from — oo to 0. Substituting y — > —y and u — > —u, and changing the limits of integration gives 



VDetM f°° , f°° , ,. N ^ 1 



pf{x,v) = - — --^ I dy I duf(-y,-u)^2^(-a x yxy~a v yvy~a xu xu-a vu vu) r 

x exp[-(a Ix i: 2 + Oj;„ot + a„„u 2 )] . (A12) 



n=0 



Due to the symmetry of the system, f(y,u) = f(—y,—u), so changing a.y — > —ay (i = j = y, u) in the matrix 
method gives alternating persistence. 



APPENDIX B: SERIES FOR p FOR RANDOM WALK AND RANDOM ACCELERATION. 

The coefficients for the random walk are: 

p(a) = 0.500000 + 0.318310 a 1 + 0.115668 a 2 + 0.021446 a 3 + 0.015651 a 4 + 0.015762 a 5 + 0.000050 a 6 - 0.003320 a 7 + 
0.007597 a 8 + 0.007587 a 9 - 0.004372 a 10 - 0.005964 a 11 + 0.004840 a 12 + 0.007913 a 13 - 0.001290 a 14 - 0.006184 a 15 - 
0.000369 a 16 + 0.004205 a 17 + 0.001467 a 18 - 0.000757 a 19 + 0.000989 a 20 + 0.000385 a 2 l - 0.003272 a 22 - 0.00 2 4 08 a 23 + 
0.003354 a 24 + 0.004972 a 25 - 0.00 72 a 26 - 0.003737 a 27 - 0.001786 a 28 + 0.00 6 50 a 29 + 0.000270 a 30 + 0.00 04 1 5 a 31 + 
0.002188 a 32 + 0.001629 a 33 - 0.001600 a 34 - 0.002448 a 35 + 0.000213 a 36 + 0.00 1451 a 37 - 0.000501 a 38 - 0.00 1 3 72 a 39 + 
0.000740 a 40 + 0.002125 a 41 + 0.000743 a 42 - 0.000602 a 43 - 0.000379 a 44 - 0.000732 a 45 - 0.001956 a 46 - 0.001143 a 47 + 
0.001780 a 48 + 0.002824 a 49 + 0(a 50 ). 

The coefficients for the random acceleration are: 

p(a) = 0.500000 + 0.477464 a + 0.021519 a 2 + 0.000314 a 3 + 0.035886 a 4 - 0.063298 a 5 + 0.029548 a 6 + 0.032884 a 7 - 
0.061472 a 8 + 0.020790 a 9 + 0.030977 a 10 - 0.039237 a 11 + 0.016107 a 12 + 0.012402 a 13 - 0.035066 a 14 + 0.021396 a 15 + 
0.020271 a 16 - 0.027422 a 17 + 0.002006 a 18 + 0.005649 a 19 - 0.009557 a 20 + 0.0 1 83 7 a 21 + 0.002099 a 22 - 0.02 3 9 60 a 23 + 
0.007071 a 24 + 0.0 1 24 05 a 25 + 0(a 26 ). 

The coefficients in both series have been truncated to six decimal places for brevity. In the plots presented in the 
paper, sufficient precision has been retained in the coefficients to ensure the accuracy of the plots and of any quoted 
values for p and 9. 
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FIG. 2. Plot of 0(a) vs a (a = e~ MAT ), where = ln[p(o)]/(21n |o|) and a < corresponds to alternating persistence. For 
a — ► 1, the series has not yet converged as the l/ln(a) term amplifies the small numerical error in p(a), whilst for a — » — 1, 
~ l/ln(|a|). 
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FIG. 3. Plot of eigenfunctions q(x) for a = 0.5 (lower curve) and a — 2.0 (upper curve, abscissa=5a;) 
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FIG. 5. Plot of /9(a), a = e~ aAT , a = 1/2 in all cases, with a < denoting alternating persistence. For a > the curves are, 
from the bottom: random walk ((3 = oo) to 0(a 49 ), and constrained Pades (3 = 24, (3 = 6, /? = 3, /3 = 2, /3 = 3/2, /3 = 1 to 
0(a 25 ). For a < the ordering is reversed, the random walk is 0(a 49 ), the random accelerations are 0(a 19 ). The Paded curves 
shown are the averages of the diagonal Pade (for odd orders) and the four nearest off-diagonal Pades, although these do not 
differ visibly (except for the plots of 8(a) for a larger than about 0.99, see figure |^). Only Pades which do not have spurious 
poles on the real axis in the range —0.25 < a < 1.25 are considered. Note that for (3 finite, p(— 1) ~ 0.108 while for the random 
walk p(—l) ~ 0.280 .... As (3 — * oo, the turn down occurs closer to a = —1 and is sharper, and so for j3 = 6, 24, the Pade can 
no longer 'predict' the curves and in these 2 cases only the raw series are plotted. 
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FIG. 6. Plot of 9(a) (= §ln(p)/ln(a)), a = e~ aAT , a = 1 /2 in all cases, where p is the same as for figure ^. a < denotes 
alternating persistence. For a > the curves are, from the top: random walk (/3 — oo), f3 = 24, f3 = 6, /3 = 3, P = 2, /3 = 3/2, 
/9 = 1. For a < 0, the order is reversed. Note that the random walk series has not converged for a — » 1, but is 1/2 in this limit. 
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FIG. 7. The Eigenfunction for a = 1/2, (3 = 3/2 and a = 0.1. 
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FIG. 8. Plots of the underdamped noisy SHO persistence exponent against a = e _AT ^ 2 , calculated using the correlator 
expansion for the cases 7 = 1/2 with lo — 2,4, 8, 16 from the lower right respectively. Also plotted are the discrete random walk 
(lowest dashed curve) and alternating discrete random walk (upper dashed curve) which are equal to the oscillating exponents 
for lo&T = 2m.iv and wAT = (2m + 1)tt respectively, m integer. Note that the series have not converged for a — » 1. 
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